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ABSTRACT 

We present non-Gaussianity testing on recently-released derived maps from the first-year WMAP data 
by Tegmark, de Oliveria-Costa and Hamilton. Our test is based on a phase mapping technique which has 
the advantage of testing non-Gaussianity at separate multipole bands. We show that their foreground- 
cleaned map is against the random-phase hypothesis at all 4 multipole bands centered around I = 150, 
290, 400 and 500. Their Wiener-filtered map, on the other hand, is Gaussian for £ < 250, and marginally 
Gaussian for 224 < £ < 350. However, we see the evidence of non-Gaussianity for £ > 350 as we detect 
certain degrees of phase coupling, hence against the random-phase hypothesis. Our phase mapping 
technique is particularly useful for testing the accuracy of component separation methods. 
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1. INTRODUCTION 

With the first-year data release of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP ), it has been pro- 
claimed that we have entered the era of "precision cos- 
mology" . The temperature fluctuations of the cosmic mi- 
crowave background (CMB) radiation are believed to be 
the imprint of primordial density fluctuations in the early 
Universe which give rise to the large-scale structures we see 
today. Hence the data enable us to test the statistical char- 
acter of the primordial fluctuations, making subsequent 
inferences on the topology and content of the Universe. 

Although the WMAP team (Komatsu et al. 2003) 
claims that the signal is Gaussian with 95% confidence, 
the internal linear combination map released by the 
WMAP team is not up for CMB studies due to "complex 
noise properties" 5 . Another group led by M. Tegmark 
has performed an independent foreground cleaning from 
the first-year WMAP data and made public their whole- 
sky CMB maps. Their foreground-cleaned map (hereafter 
FCM) and the Wiener-filtered map (WFM) are available 
online (Tegmark 2003). 

The FCM by the authors' definition is such that the 
foreground contamination is removed as much as possi- 
ble. As foregrounds are rather non-Gaussian, any residual 
after cleaning would manifest itself in the phase configu- 
ration. In this Letter we display the phases of the FCM 
and the WFM with color coding and implement our phase- 
mapping technique to test quantitatively the Gaussianity 
of both maps, based on the random-phase hypothesis of 
homogeneous and isotropic Gaussian random fields. Our 
phase mapping technique can play a crucial role as a qual- 
itative criterion for component separation similar to the 
field of image reconstruction. 



2. GAUSSIAN RANDOM FIELDS AND THE RANDOM PHASE 
HYPOTHESIS 

The statistical characterization of temperature fluctua- 
tion of CMB radiation on a sphere can be expressed as a 
sum over spherical harmonics: 



AT(0, <p) 
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where at m = \ai m \ exp(i<f>t m ). Homogeneous and isotropic 
Gaussian random fields (GRFs), as a result of the sim- 
plest inflation paradigm, possess Fourier modes whose 
real and imaginary parts are independently distributed. 
In other words, they have phases 4>i m that are indepen- 
dently distributed and uniformly random on the interval 
[0,2tt] (BBKS 1986; Bond & Efstathiou 1987). Thus the 
spatial variations should constitute a statistically homoge- 
neous and isotropic GRF (BBKS 1986) whose statistical 
properties are completely specified by its angular power 
spectrum Ci, 

(a em a*, m , ) = Ce 6 u >5 mm i . (2) 
The strict definition of a homogeneous and isotropic 
GRF requires that the amplitudes are Rayleigh distributed 
and the phases are random (Watts & Coles 2003). At the 
same time, the Central Limit Theorem guarantees that a 
superposition of a large number of Fourier modes with ran- 
dom phases will be Gaussian. Therefore the random-phase 
hypothesis on its own serves as a definition of Gaussianity 
(BBKS 1986). 

3. COLOR-CODED PHASE MAP OF THE DERIVED 
WMAP MAPS 



Tegmark, de Oliveira-Costa & Hamilton (2003) 
(TDH03) perform an independent foreground analysis 
from the WMAP data and provide a FCM and WFM. 
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We first use a visual display of phases by colors to show 
phase associations (Coles & Chiang 2000). In color im- 
age display devices, each pixel represents the intensity and 
color at that position in the image. Two color schemes 
are usually used for the quantitative specification of color, 
namely the Red-Green-Blue (RGB) and Hue-Saturation- 
Brightness (HSB) color schemes. Hue is the term used to 
distinguish between different basic colors (blue, yellow, red 
and so on). Saturation refers to the purity of the color, 
defined by how much white is mixed with it. Brightness 
indicates the overall intensity of the pixel on a grey scale. 
The HSB color model is particularly useful because of the 
properties of the 'hue' parameter, which is defined as a 
circular variable. Therefore we are mapping phases from 
to 2-7T to the hue circle. 

We have used the Healpix 6 package to produce ai m . 
In Fig. 1 we show the color-coded phase gradient Dg = 
4>i+i,m — <j>t, m f° r the FCM and WFM. The vertical axis 
is the multipole I up to £ — 600 and the horizontal the 
m axis where m < £. Due to the relation ag m — a\ _ m , 
only modes from non-negative m are shown. Although 
the phase gradient (from neighboring modes) is the most 
primitive way of qualitatively checking phase correlations, 
the apparent presence of stripes shown in the FCM indi- 
cates strong coupling between modes of neighboring £ of 
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Fig. 1. — Color-coded phase gradient Di for the FCM (upper 
left triangle) and the WFM (bottom right). The vertical axis is 
the I up to I = 600 and horizontal the m axis. Due to the relation 
a l,m = a i _ m i we only show modes from non-negative m. Although 
the phase gradient (from neighboring modes) is the most primitive, 
the stripes shown from the FCM indicate strong phase correlation 
between modes of neighboring I of the same m. 

4. PHASE MAPPING AND THE MEAN CHI-SQUARE 
STATISTIC OF THE DERIVED MAPS 

To test the Gaussianity of the FCM and the WFM based 
on the random-phase hypothesis, we apply a phase map- 
ping technique (Chiang, Coles & Naselsky 2002; Chiang, 
Naselsky & Coles 2002) to quantify the degree of 'ran- 

6 http : //www . eso . org/science/healpix/ 



domness' of the phases (i.e. Gaussian). The return map 
of phases is a bounded square in which all phase pairs 
of fixed separation (Am, A£) are mapped as points (see 
Fig. 2). For example, one single return map for phase 
pairs with separation (Am, A£ ) = (0, 1) contains points 
with (x,y) coordinate {<pi. m , 4>i+i,m), i-e. all phase pairs 
from modes that are separated by A£ = 1. If the phases 
are random, we expect to have an ensemble of return maps 
of all possible separations, each of which should be a scat- 
ter plot. As such we are testing the 'randomness' on the 
most strict terms. After mapping phase pairs on to a re- 
turn map, we can apply a mean \ 2 statistic on the return 
map, which is defined as 

— 1 \ - [pis j) -f\ ,„s 

x= m22 p ' W 

where M is the number of pixels on the return map, p 
is the mean value for each pixel on the discretized return 
map. Chiang, Naselsky & Coles (2002) have shown that 
for a homogeneous and isotropic GRF, return mapping of 
phases results in an ensemble of return maps, each with 
a Poisson distribution. The expectation value of the \ 2 
from such ensembles of Poisson-distributed maps is 

(Bp = ^ (4) 

where R is the scale of smoothing from a 2D Gaussian 
convolution in order to probe the spatial structure. The 
X 2 P will have a statistical spreading around (x 2 )p with a 
dispersion Ep where 
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Fig. 2. — An example of a return map for (Am, A!) = (0,4) of 
phases (pe,m °f the FCM where 41 < £ < 250. The \ 2 of this return 
map is 0.0332 when it is discretized into 128 2 pixels with smoothing 
scale R = 2. 

Figure 3 shows the histograms of the x 2 statistics from 
the ensemble of the return maps of the FCM and the WFM 
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for 4 multipole bands. One of the advantages of phase 
mapping technique is that we are able to check Gaussian- 
ity in different multipole bands, in particular those cor- 
responding to foreground contamination and noise. Here 
we present the x 2 statistic at 4 bands centered around 
£ ~ 150, 290, 400 and 500: 41 < £ < 250 (roughly the 
first Doppler peak), 224 < £ < 350, 350 < £ < 450 and 
463 < £ < 550. The solid dark and dotted gray curves are 
the WFM and FCM, respectively. In each panel the verti- 
cal line denotes the expectation value (x 2 )p = (47ri? 2 ) -1 . 
The curves from the FCM are obviously skewed, and 
hence are manifestations of phase correlations (i.e. non- 
Gaussian). 
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Fig. 3. — The histograms of x 2 statistic for the FCM (the dot- 
ted gray curves) and the WFM (the solid dark curves) at different 
multipole ranges I. One of the advantages of the phase mapping 
technique is that it enables us to check non-Gaussianity for dif- 
ferent multipole ranges. The upper horizontal axis is annotated in 
terms of the theoretical dispersion Sp of GRFs with origin set at the 

expectation value (x 2 )p = (iTrR 2 )^ 1 (vertical line in each panel). 
The smoothing scale on the M = 128 2 discretized return map is 
R = 2. 

In Fig. 4 we display the gross behavior of the distribu- 
tion curves in terms of the arithmetic mean (% 2 ) and the 
dispersion S from the mean chi-square statistic. The top 
panel is from the FCM and the bottom WFM. The con- 
tours mark 68% (solid curve) and 95% (dotted curve) CL 
regions from 2000 realizations of GRFs. The symbols cor- 
responds to 4 multipole bands centered at £ ~ 150, 290, 
400 and 500. Note that the contour region in the bot- 
tom panel corresponds to a small section in the top panel. 
The phases of the 4 multipole bands from the FCM are all 
strongly correlated, so they are far away from the 95% CL 
region. The WFM, however, shows that phases below the 
first Doppler peak are random, with the other 3 multipole 
bands around the edge of 68% CL region. 

We see evidence of non-Gaussianity, however, in the 
WFM of the following two bands centered £ ~ 400 and 500. 
In the lower two panels of Fig. 3 there are points appearing 



at the tails above 6£ P . On the other hand, among the 2000 
realizations we simulate for GRFs, no mapping of phases 
reaches x 2 value over 6E P , setting the probability below 
0.05% for a GRF to have such mapping. Phase mapping 
from the separation (Am,A£) = (0,2) produces x 2 value 
at 7.3Sp for the multipole band centered £ ~ 400, also at 
6.5S P at (Am, A£) = (1,2). For the band £ ~ 500, 7.6S P 
appears at (Am, A^) = (2,2). These phase couplings are 
clear signs against the random-phase hypothesis, therefore 
a manifestation of non-Gaussianity. 
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Fig. 4. — The mean chi-square statistic from the FCM (top) and 
WFM (bottom) against 2000 realizations of GRFs, which is dis- 
played in terms of the arithmetic mean (\ 2 ) an d the dispersion E 
of their distribution curves. The contours mark 68% (solid curve) 
and 95% (dotted curve) CL regions from 2000 realizations of GRFs. 
Although 68% and 95% denotes 1 and 2-cr deviation in Gaussian 
statistics, the distribution is not Gaussian but rather chi-square. 
The cross (x), triangle (A), square (□) and star (*) symbols denote 
X 2 statistic from multipole ranges centered I ~ 150, 290, 400 and 
500, respectively. Note that the contour region in the bottom panel 
corresponds to a small section in the top panel. 

We plot in Fig. 5 the CMB temperature map from 
only two multipoles £ — 350 and 352 (of all m's) of the 
FCM and WFM. The choice of these specific multipoles 
of A£ — 2 from our previous calculation is to demonstrate 
non-Gaussian signals the correlated phases will produce in 
the map. The structures at tp ~ and it in the FCM, the 
residual signal after foreground cleaning, disappear after 
Wiener filtering. 
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Fig. 5. — The CMB temperature from two multipoles £ = 350 
plus 352 of the FCM (top) and the WFM (bottom). These two 
multipole modes are chosen because of the pronounced coupling be- 
tween modes A£ = 2 of all m's. The structures at ip ~ and tt 
shown in the FCM disappear after Wiener filtering, from which it 
is marginally Gaussian at these two multipoles. 

5. DISCUSSIONS 

In this Letter we have tested non-Gaussianity of two 
maps: the Foreground-cleaned map and the Wiener- 
filtered map, which are processed by TDH03 from the 



WMAP data. Based on the random-phase hypothesis, we 
use a phase mapping technique to yield a statistic that has 
detected considerable non-Gaussian signals for both maps 
at most multipole bands. Our phase mapping technique 
is particularly useful in separating non-Gaussian contribu- 
tions from different sources when various contaminations 
are present at different £ ranges. A multipole band which 
is considerably non- Gaussian could have an insignificant 
non-Gaussian contribution in the whole map and still pro- 
duce an over-all Gaussian realization within a certain con- 
fidence level. As the uncertainties in foreground cleaning 
propagate through the data-processing pipelines to the ac- 
curacy of the angular power spectrum, it is therefore nec- 
essary to have effective methods in component separation. 
We believe that our phase mapping technique is a use- 
ful criterion to be incorporated into such methods. Our 
statistic based on phase mapping also holds great advan- 
tage when it comes to the issue of creating many whole- 
sky Gaussian realizations for Gaussian statistics. As our 
null hypothesis is that phases are random, we only need 
to put random phases (with Gaussian instrumental noise 
being automatically included) for each harmonic mode, 
which is easily done without any limit on the highest har- 
monic number I from any pixelization scheme. It is worth 
mentioning that the upcoming Planck mission will have 
higher sensitivity and resolution, hence every step of data 
processing will be crucial in reaching such precision. 
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